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The relativistic fermion-antifermion bound state vector potential of constraint theory is 
calculated, in perturbation theory, by means of the Lippmann— Schwinger type equation 
that relates it to the scattering amplitude. Leading contributions of n-photon exchange 
diagrams are calculated in an approximation scheme that adapts eikonal approximation to 
the bound state problem. They produce terms proportional, in three-dimensional x-space, 
to (a/r) n . The series of leading contributions is summed. 
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We report in this note results on calculations, in the framework of the manifestly 
covariant formalism of constraint theory |2|, |3], of the relativistic two-body bound state 
vector potential from summation of Feynman diagrams, considered in their leading order 
of the infra-red counting rules of QED. 

Perturbative calculations of the three-dimensional potential from the Bethe— Salpeter 
kernel or from the scattering amplitude have been successful in the past mainly in the 
Coulomb gauge for the exchanged photon propagator Q. Covariant gauges ||, or covari- 
ant propagators for scalar exchanges , produce spurious infra-red singularities that are 
only cancelled by contributions of higher order diagrams. This defect seems to be related 
to the noncovariant nature of the three-dimensional equation around which iteration of 
the Bethe— Salpeter equation is accomplished. It turns out, as we shall see below, that in 
the manifestly covariant formalism of constraint theory these complications are absent: in 
each formal order of perturbation theory the leading term of the potential due to ladder 
and crossed ladder diagrams is free of spurious singularities. A similar phenomenon also 
occurs in the variant of the quasipotential approach developed by Todorov 0. In partic- 
ular, it has been shown ||, for one spin- 1/2 and one spin-0 particle systems, that in the 
two-photon exchange diagrams, calculated in the Feynman gauge, the spurious infra-red 
singularities cancel out. 

We concentrate in the following on the case of bound states composed of one fermion 
with mass m\ and one antifermion with mass 7712, mutually interacting by means of vector 
(massless) photons. 

We use standard notations for the total and relative variables: P = pi + P2, P = 
(pi — P2V2, x = x\ — X2- For states that are eigenstates of the total momentum P we 
define transverse and longitudinal decompositions of four-vectors with respect to P and 
denote them with indices T and L, respectively. A detailed presentation of the constraint 
theory wave equations, of their properties, as well as of the way they can be reduced to 
a single Pauli— Schrodinger type equation can be found in Ref. ||. The wave equations 
determine, among others, the cm. relative time evolution of the system and allow the 
elimination of the cm. relative energy. These properties are consequences of the following 
constraint equation: 

C(p) = {pi - p\) - (ml - m 2 2 ) = 2P L p L - (m? - m\) « . (1) 
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When constraint C is used, the individual Klein— Gordon operators become equal: 



H = (p\-m\) 
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The potential V that appears in the wave equations must be compatible with constraint 
C [Eq. (HDJ. This implies that it should not depend on the relative longitudinal coordinate 
xl- V = V(x T , Pl,p t , 71, 72) • [71 and 72 are the Dirac matrices for particles 1 and 2, 
respectively. In general, V is an integral operator in x T .} Thus, the internal dynamics of 
the system is three-dimensional, apart from the spin degrees of freedom. 

The wave equations of constraint theory can be connected to the Bethe— Salpeter 
equation, by iterating the latter around the constraint hypersurface ([I]) [IU|. This con- 
nection determines the relationship of the potential V with the off-mass shell scattering 
amplitude, through a Lippmann— Schwinger type equation: 



V = T + Vg T , 
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where T is the off-mass shell scattering amplitude in which the external momenta are 
submitted to the constraint ([I]); g is defined as: 
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Si{pi)S 2 (-p 2 )H 



C(p) 



(4) 

to the con- 



where 5*i and S 2 are the propagators of the fermions (submitted in Eq. 
straint ([[])) and H is defined in Eq. (0). 

Equation (||) is the basis for the calculation of the potential in terms of Feynman 
diagrams. We limit ourselves to the evaluation of the ladder and crossed ladder diagrams 
(in their leading order), neglecting the contributions of radiative corrections. Iteration of 
Eq. (Rt) yields for the potential the expansion: 



p=0 



(5) 



71=1 



where V^ n ' is that part of the potential that results from n exchanged photons. We observe 
that the perturbation series of V contains, in addition to the usual Feynman diagrams, 
other types of diagram arising from the presence of the constraint factor g ; we shall call 
these diagrams "constraint diagrams" ; they are obtained from the usual box-ladder type 
diagrams by the replacement of (at least) one pair of fermion and antifermion propagators 
by the constraint factor 2m5(C)go. The role of these diagrams will be to cancel the 



spurious infra-red singularities coming from the amplitude. The diagrams corresponding 
to two photon exchanges are represented in Fig. 1. [The constraint diagram is denoted 
by a cross.] 

When constraint ([!]) is imposed on the external particles, one obtains: = Pil—Pil = 
0, q 2 = (pi — p[) 2 = q T2 . The Feynman diagrams are calculated with external particles 
considerd off the mass shell, with their longitudinal momenta fixed by the bound state 
mass Pl [the binding energy is of order a 2 }, through Eq. ([]]), while the transverse momenta 
p T and p' T have, according to the infra-red counting rules of QED, orders of magnitude 
of a. The photon is taken massless. 

The potential that results from one photon exchange is in the Feynman gauge, in 
three-dimensional x-space: 

y {1) = 2^7i.7> , a = r = (6) 

[The matrices 72, corresponding to the antifermion, act on the 4x4 matrix wave function 
from the right.] 

The vector potential can be divided into two parts: the timelike component, propor- 
tional to 71L72L, and the spacelike component, proportional to 71^72^- Because of the 
latter matrices, which are of order a 2 , the spacelike component, at each formal order of 
perturbation theory, will be damped with respect to the timelike component by an 0(a 2 ) 
factor. In the present work, where we are calculating only leading terms of higher order 
diagrams, we can evaluate only the timelike component of the vector potential. However, 
if one assumes that in the Feynman gauge a simple structure of the 7-matrices like that 
of Eq. @ survives in higher orders, then the knowledge of the timelike component will 
still allow one to reconstitute the whole potential in this gauge. 

The potential is calculated by evaluating the diagrams of Fig. 1. The details of 
the calculations will be presented elsewhere; we also refer the reader to Ref. ||, where 
analogous calculations are done and similar properties to those obtained here found. The 
various integrals can be calculated to the desired precision by dividing the integration 
domains into several intervals, in which appropriate approximations can be used. 

To understand the various cancellation mechanisms, we explicitly write the product 
of the two fermion propagators occurring in the box diagram of Fig. 1, where k is the 
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four-momentum carried by one of the photons: 
Si{ P i-k)S 2 {-(p 2 + k)) 



(pi — k) 2 — m\ + ie (p 2 + k) 2 — m\ + it 
(li.px + mi)(-7 2 .p 2 + m 2 ) - (ji-Pi + mi)j 2 .k - (-j 2 .p 2 + m 2 )ji.k + j x .kj 2 .k 



(7) 



The first term in the brackets yields a contribution to the potential (in x-space) com- 
posed of two terms, of orders a 2 In a -1 and a 3 In a -1 , respectively The analogous con- 
tribution of the crossed diagram is of order a 3 In a" 1 . The corresponding contribution 
of the constraint diagram is of order a 2 In a -1 and cancels the similar term of the box 
diagram. The remaining term of the box diagram cancels the 0(a 3 lna _1 ) term of the 
crossed diagram and yields a term of order a 3 , which in turn can be removed by a finite 
multiplicative renormalization of the constraint factor (^) [we use the prescription that 
the 1 / r-terms of the potentials should arise only from the one-photon exchange diagram] . 
The sum of the above terms is of order a 5 In a -1 , which, compared to the leading 0(a) 
effect we are estimating, can be neglected. 

The second and third terms in the brackets of Eq. (0) and the similar terms of the 
crossed and constraint diagrams yield 0(a 4 lna _1 ) and 0(a 4 ) terms. The 0(a A \aa~ 1 ) 
terms are mutually cancelled and one remains with the 0(a A ) terms. The latter actu- 
ally arise from the longitudinal component ki contributions of the photon momentum k 
multiplying the 7-matrices. 

The last term in the brackets in Eq. (|7|) and the similar ones of the other diagrams 
are individually negligible in front of the 0(a 4 ) terms. 

For the evaluation of the 0(a A ) terms, one can replace at leading order the matrices 
7ii and 7 2 l by their eigenvalues, +1 and —1, respectively, with respect to the dominant 
component of the wave function, as well as mi by p 1L and m 2 by p 2L . One finds: 

y (2) = -^27lL72L • (8) 
L 

It can be shown || that, when this contribution is incorporated in the wave equations, 
the latter yield to order 1/c 2 , up to wave function transformations, the Breit Hamiltonian, 
which is known to produce the correct 0(a A ) effects. This means that no other 0(a 4 ) 
terms should arise from higher order diagrams. 
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It is natural to expect that the cancellation mechanism that occurs at the level of 
two-photon exchange diagrams and yields at leading order a local potential generalizes 
to higher orders. It does not seem possible to verify this property in a rigorous way, due 
to the complexity of the structure of n-photon diagrams. However, within a plausible 
approximation scheme, that takes into account the properties of the bound state we are 
considering, it is possible to show the above property. 

To device the appropriate approximation scheme, we make the following observations: 
i) In result @ the factor 1/r 2 represents in momentum space the three-dimensional con- 
volution (up to multiplicative coefficients) of the two photon propagators, in which the 
longitudinal component k\ has been neglected; this means that the entire q 2 dependence 
of the potential (at leading order) comes from the photons and not from the fermions. ii) 
is independent of the off-mass shell condition imposed on the fermions (i.e., does not 
depend on (p\ — ml) or [j>\ — m|)). iii) In the product of fermion propagators [Eq. ([?])], 
only the linear longitudinal component terms fez, of the numerator have contributed to 
V^ 2 \ These observations suggest the following approximation of the fermion propagators: 

Si{pi - fei) ^ — — : — — [{llLPlL + mi) - 7li*li] > 
—zpi.Kx + ie 

S 2 (-(P2 + h)) - ~ ; — — [(-72LP2L + m 2 ) - liikiiX ■ (9) 

In the numerators, we have neglected the transverse momenta p T and kj (but not in 
the denominators); in products of two fermion propagators we neglect in the numerator 
quadratic terms k\ (of the same fez); we also neglect the total momentum transfer q. In 
photon propagators, we neglect k\ dependences when they appear in nonleading terms. 
Approximation (|9|) may be considered as a variant of the eikonal approximation |]TT , |12 



adapted to the bound state problem. Finally, because of the on-mass shell treatment of 
the external fermions in approximation (|9D, the photon should be given a small mass to 
prevent infra-red divergences at intermediate stages. 

It can be checked by direct calculation that the above approximations globally produce 
the correct result (|8]). 

We now apply approximation (^) to n-photon exchange diagrams (n > 2). The main 
formula we use is a generalization of an identity already used in the eikonal approximation 
(cf. Ref. Jl2|, Appendix). Let (c 1; c 2 , • • • , c n ) be a set of n numbers; it can be divided 
into n independent subsets of (n — 1) numbers, where in each subset one of the q's 
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n) is missing. We have the identity: 



i=l perm 



EE (c' 1 )-\c' 1 + c' 2 r 1 ...(c' 1 + c' 2 + --- + c' n _ 1 )- 1 = X> / II<* , (10) 



n \ i I n 



u=l / ' W=l 



where (c' 1; c^,, • • • , c^_ x ) is a permutation of the subset of (n — 1) numbers (ci, c 2 , • ■ ■ , Cj_i, 
Cj+i ; * ' ' j Cn) • 

Let us consider an n-photon exchange diagram (Fig. 2). We denote by k\, ■ ■ ■ , k n the 
momenta carried by the photons. To take into account momentum conservation, a factor 
(27r) 4 5 4 (X]™ = i hi ~ o) mus t be included in the corresponding integral. Since q L = 0, we 
have YJl=\ kiL = 0; but according to the approximations (H) we may also neglect q T in the 
fermion propagators; hence, we have there the approximation 

n 

£fe = 0. (11) 

8=1 

The total number of the diagrams considered above is n\. In a given diagram we have 
(n — 1) propagators of fermion 1 and (n — 1) propagators of fermion 2; only (n — 1) 
photon momenta appear on each of the fermion lines (but not the same in general). We 
first consider the terms that do not contain any kn (i = l,---,n) in the numerator 
[Eq. (§)]. The corresponding n! diagrams can be divided into n sets, where in each set 
we have (n — 1)! permutations of fermion 2 propagators (without the kn factors in the 
numerator) containing the same set of (n — 1) photon momenta. With approximations 
(H), the sum of all these propagators has the structure of the left-hand side of Eq. (|T0|), 
with Cj = (2p 2 .ki + ie) [the terms i{— 72LP2L + ^2) are factorized]; hence, we obtain the 
factor 



n-l 



'2,n — 



Y[(2p 2 .kj + if) 



-2t7T) n - 1 5(2p 2 .k 1 )5(2p 2 .(k 1 + k 2 )) ■ ■■5(2p 2 . ^ h) , (12) 



i=i 



3=1 



the second equality resulting from the use of Eq. (JTTJ) and several algebraic operations (for 
more details cf. Ref. pp. 116-117 and also Ref. ||11|| ). The 5-functions can then be 
used, upon integrations on the fc^'s, in the fermion 1 propagators to yield a global factor 
{2PL)~^ n ~ l \^iLpiL + mi) n ^ 1 {— r )2LP2L + rn2) n ~~ l ■ Denoting by the total contribution of 
the above diagrams (without the coupling constants and other multiplicative coefficients), 
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we obtain: 

/ — i \ n ~ 1 

J in) = {^-jrj [(IilPil + ■m\){-l2LViL + m 2 )} n ~ 1 (7il72l)™ 

1 1 l 

X (-2p T .kf + ie) (-2p T .(kf + kl) +iej"' (-2p T . Eg kj + le) ' 

(13) 

where fi is a small mass given to the photon, to prevent infra-red divergence. Also, the 
5-functions fll2f) yield for the kf L -terms orders of magnitude of 0(a 4 ), which accounts for 
their omission in front of 0(a 2 )-terms in approximations (^) and in photon propagators. 

On the other hand, we should also take into account the contributions of constraint 
diagrams that are associated with the diagrams considered above. A typical diagram 
where the constraint factor go [Eq. @] appears twice is shown in Fig. 3. The analysis 
of these diagrams is very similar to that already done in the infinite mass limit (cf. Ref. 
P, Appendix C). Each constraint diagram contributes the same quantity above 
up to a sign factor equal to (— l) p if g occurs p times in expansion (|5|). Taking into 
account the combinatorial factors one finds that the sum of constraint diagrams cancels 
the contribution of ordinary diagrams. 

This result can be generalized to the contributions of those parts of the above diagrams 
where some of the kn's (but not all) appear in the numerator. These parts have also 
constraint diagram counterparts and are cancelled by them. 

Therefore, the only surviving parts of the above diagrams are those where products 
of (n — 1) independent combinations of the fc^'s appear in the numerator; these parts 
do not have constraint diagram counterparts, since there kn = (i = 1, • • • ,n). After 
the appearance of the (^-functions from the permutational procedure [Eq. (|12|)1, realized, 
according to the origin of the /c^l's, partly on line 1 and partly on line 2, the fc^'s can- 
cel by integration, up to multiplicative coefficients, the remaining fermion propagators, 
and one obtains a three-dimensional convolution of n photon propagators, yielding in x- 
space a potential proportional to (a/r) n . The analysis is much simplified by using these 
cancellations (with the appropriate coefficients) prior to the permutational procedure. 

We next turn to the calculation of the combinatorial factor associated with the above 
procedure. In the product of the (n — 1) pairs of fermion and antifermion propagators of 
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a given diagram we can select p fc^'s on line 1 and the complementary (n — l—p) fciz's on 
line 2. The choice of p k^s on line 1 can be made in ~p~ ^ different ways; however no 
freedom is left for the choice of the set of (n — 1 —p) complementary fc^'s on line 2. During 
the permutational procedure, we need (p+l)\ permutations to obtain p ^-functions on line 
2 from the propagators not having fc^'s in their numerator; similarly, we need (n — p)\ 
permutations to produce (n — l—p) 5-functions on line 1. After the approximations 
(IilPil + wii) — 2piL and (— J2LP2L + 1^2) — %P2L are made, the integrations on the /cjl's 
produce the global factor ((— 2m) / (2P L )) n ^ 1 . 

Taking into account the total number of the diagrams with n exchanged photons, 
which is n!, we obtain for the combinatorial factor: 



c - t C p ') ( P+1) t- P) ! = F ( -<" - -« 2; ^ = jJ+Tjw. • (14) 



p=0 



where F(a, b; c; z) is the hypergeometric function [T^|, the value of which for z — 1 is 
known in terms of T-functions. 

The corresponding potential then becomes: 

V (n) = —(-ie) 2n i 2{n - l) (-i) n ( — V"' (2n)! ( llLl2 L) n 
2P L K ' 1 ; \2P L ) n + 1 )\n\ UlL72L) 



yielding in x-space: 



y<»> = 7ll72i . (i 6 ) 



(n+l)!n! V 2P L r 



[We have replaced (7il72l)™ _1 by ( — l) n_1 .] 

This expression is also valid for the one-photon exchange case [Eq. (|6])], for which it 
provides the exact result, concerning the timelike component. 

The total timelike potential becomes: 

v = £(-ir (2n)! f— V* = -mm ( . (17) 

A parametrization that allows external field interpretation of the potentials and en- 
sures positivity of the norm for energy independent potentials was proposed by Crater 
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and Van Alstine and used in Ref. ||; it is: V = tanhV. Potential V above fits this 
structure. Writing 

V = tanh( 7lL72L y 2 ) , (18) 

we identify 

* = M 1 + !;) • < 19 > 

Furthermore, if we adopt, for the vector potential in the Feynman gauge, the structure 
(HD of the 7- matrices, obtained in lowest order, by the substitution 71x721,— ►Ti-Ta, we can 
reconstitute the whole potential in this gauge: 

V = tanh( 7l .7 2 V 2 ) , (20) 

with V 2 given by Eq. ([191) . Expression (|2~0D is nothing but the fermionic generalization 
of the electromagnetic potential proposed by Todorov in the quasipotential approach for 
two spin-0 particle systems, on the basis of minimal substitution rules M , and later used 



in spectroscopic calculations [g, |15| . It should be stressed that although potential V 
above is bounded in the limit r — > 0, the corresponding effective potential of the final 
Pauli— Schrodinger type equation of the reduced wave function |J has the usual 1/r and 
1/r 2 singularities. Furthermore, in the infinite mass limit for one of the particles, it is only 
the one-photon exchange contribution that survives in the wave equations and the latter 
reduce to the Dirac equation of the finite mass particle in the presence of the Coulomb 
potential. 

The major result of the present work is that the sum of the leading terms of multipho- 
ton exchange diagrams provides a local expression for the potential in three-dimensional 
x-space. This feature considerably simplifies the analysis and resolution of the corre- 
sponding wave equations, without loss of the relativistic invariance of the theory, neither 
of the effects of multiphoton exchanges. These effects become more sizable in strong cou- 
pling problems, such as evaluation of the short-distance part of the interquark potential 
in QCD, which is very similar to the QED potential, with values of the coupling constant 
of the order of 0.5, or probe of the strong coupling regime in QED. Formal generaliza- 
tions to treat with the present summation method effective interactions, in particular the 
confining ones, are possible as well. The cases of spin-0 particles and of scalar interaction 
and a more detailed account of this work will be presented elsewhere. 
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Figures 



Fig. 1. Two-photon exchange diagrams; the "constraint diagram" is denoted by a cross. 
Fig. 2. A multiphoton-exchange diagram. 

Fig. 3. A typical "constraint diagram" where the constraint factor go [Eqs. (f|) and (||)] 
appears twice (p = 2). 
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